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Представлена структура программного комплекса для исследования характеристик термоэмиссионных систем. Описан процесс 
расчета характеристик электрогенерирующей сборки, состоящий из трех этапов: подготовительный, расчетный и визуализации 
результатов. В качестве иллюстрации возможностей программного комплекса приведены результаты расчета аномальных ха- 
рактеристик электрогенерирующих сборок. 


Практическая реализация технологии вычисли- 
тельного эксперимента базируется на трех соста- 
вляющих: модели, алгоритмы и программы. Во 
многих случаях для проведения вычислительных 
экспериментов достаточно использовать современ- 
ные математические пакеты, в которых реализова- 
ны эффективные и проверенные временем числен- 
ные методы решения типовых задач. При модели- 


ровании сложных физических установок, к кото- 
рым, безусловно, относятся термоэмиссионные си- 
стемы, на этапе создания математического описа- 
ния производится декомпозиция взаимосвязанных 
процессов по различным критериям, что приводит 
к созданию множества иерархических моделей [ 1 ] и 
алгоритмов [2] разной степени сложности. Необхо- 
димость учета влияния на исследуемые процессы 
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конструкций и технологий изготовления физиче- 
ских установок приводит к многовариантности 
множеств иерархических моделей и, что не менее 
важно, проблеме хранения исходных данных и ре- 
зультатов вычислительных экспериментов. 

Таким образом, основными требованиями к 
программному обеспечению являются не только 
реализация эффективных численных методов в 
каждом программном модуле, но и возможность 
либо выбора из существующего множества моде- 
лей наиболее подходящей, либо оперативное соз- 
дание новой модели из библиотеки модулей для 
проведения целенаправленных вычислительных 
экспериментов, формирования, визуализации и 
хранения результатов. 

Кроме того, в связи с постоянным совершен- 
ствованием отдельных программных модулей, а 
также с целью повышения достоверности результа- 
тов программные модули, выполняющие одинако- 
вые функции, должны быть взаимозаменяемы. 

В наибольшей степени перечисленным требо- 
ваниям удовлетворяет подход, основанный на раз- 
делении функций подготовки исходных данных в 
виде банка данных, построения модели в виде це- 
почки последовательно выполняемых програм- 
мных модулей, собственно проведения расчетов, а 
также визуализации полученных данных [3]. Моде- 
лирование характеристик электрогенерирующих 
сборок (ЭГС) выполняется в три этапа [4]. На пер- 
вом этапе с помощью подготовительных модулей 
формируются файлы банка данных, в которых со- 
держится информация о геометрических параме- 
трах, электрофизических, эмиссионных и тепло- 
физических свойствах элементов ЭГС. На втором 
этапе проводится расчет, визуализация и анализ се- 
рий вольт-амперных характеристик термоэмисси- 
онного преобразователя (ВАХ ТЭП) для предпола- 
гаемых условий работы ЭГС. На третьем этапе рас- 
считываются температурные и электрические поля 
и ВАХ ЭГС. Условия расчета ЭГС (распределение 
тепловыделения, тепловая мощность, давление па- 
ра цезия в резервуаре, входная температура и рас- 
ход теплоносителя, сопротивление нагрузки и ко- 
личество расчетных точек) задаются в управляю- 
щих файлах банка данных. Для моделирования 
аномальных ситуаций предусмотрены специаль- 
ные подготовительные модули, позволяющие фор- 
мировать исходные данные, имитирующие возни- 
кновение утечек тока между электродами, через 
коллекторную изоляцию и между элементами, а 
также изменения теплофизических свойств мате- 
риалов, распухание топливно-эмиттерных блоков, 
вынос топлива в межэлектродные зазоры и т.д. 

Входные, выходные и промежуточные данные 
работы функциональных модулей хранятся в банке 
данных в виде текстовых файлов. 

Расчетные модули создаются на основе про- 
грамм (процедур и функций), которые реализуют 
математическое [1] и алгоритмическое обеспече- 
ние [2]. Вспомогательные программы также офор- 


мляются в виде модулей и объединяются в библио- 
теки: 

• стандартных численных методов; 

• программ расчета коэффициентов аппроксима- 
ции уравнений тепло- и электропроводности, 
теплофизических и электротехнических 
свойств материалов, составляющих тепловых 
потоков и т.п.; 

• чтения файлов исходных данных и записи ре- 
зультатов расчета. 

Модули визуализации позволяют просмотреть 
как исходные данные, так и результаты промежуточ- 
ных и окончательных расчетов в виде распределений 
температурных и электрических полей, ВАХ и т.д. 

Программные и функциональные модули реа- 
лизованы с помощью системы программирования 
ТигЬо Раксаі 7.0. На рис. 1 представлена структур- 
ная схема взаимодействия функциональных моду- 
лей, библиотек программ и банка данных. 



Рис. 1 . Структура программного обеспечения для модели- 
рования термоэмиссионных систем 


Для проведения вычислительных эксперимен- 
тов в качестве системного программного обеспече- 
ния используется пакет программ ВгаіпВіогт [5], 
который обладает следующими возможностями: 

• объединение функциональных модулей, реали- 
зованных в М5-008, в пакет программ с мини- 
мальными изменениями их кода; 

• добавление новых и замена существующих 
функциональных модулей; 

• запуск модулей на исполнение и автоматизация 
последовательного выполнения цепочек модулей; 

• предоставление доступа к данным и выполнение 
операций над ними: редактирование, копирова- 
ние, добавление и удаление без выхода из си- 
стемной оболочки; 

• повышение надежности работы в системной 
оболочке, в частности, аварийные остановы в 
модулях не приводят к разрушению файлов и 
оболочки. 

Базовое окно системной оболочки Вгаіп8іогт с 
загруженным пакетом прикладных программ для 
исследования процессов в термоэмиссионных си- 
стемах представлено на рис. 2. 
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Рис. 2. Базовое окно системной оболочки Втаіпзіотт 


В таблицах 1-4 приведен перечень и назначе- 
ние основных функциональных модулей, с помо- 
щью которых проводятся исследования различных 
режимов работы термоэмиссионных систем. 


Таблица 1. Состав функциональных модулей ТЭП 


Наименова- 
ние модуля 

Назначение модуля 


Подготовка: 

РСР02 

данных для зависимости давления пара цезия 
от температуры резервуара с цезием 

РСР05 

кривых Рейзора электродов 


Расчет: 

РСС01 

вакуумных и квазивакуумных ВАХ 

РСС02 

диффузионных ВАХ 

РСС03ш_р 

ВАХ по эмпирической и полуэмпирической мо- 
делям плазмы ТЭП 


Таблица 2. Состав модулей для расчета электрических про- 
цессов эгс 


Программный мо- 
дуль 

Назначение модуля 


Расчет: 

РРОІ 

удельного электрического сопротивления 
материалов 

5ХЕМА1 

коэффициентов аппроксимации уравне- 
ний электропроводности 

ѴОНЗ, 

ѴОІД4, 

ѴОИ31, 

ѴОІ_Т4І 

параметров линейных эквивалентов ТЭП 
с использованием трех или четырехмер- 
ной интерполяции 

ЕІ_ЕСАІ\І1 

распределений потенциалов, токов, на- 
пряжений и интегрального значения тока 


Верификация различных версий программных 
модулей производилась в течение ряда лет путем 
сравнения расчетных результатов с результатами 
других авторов, а также с экспериментальными дан- 
ными, полученными при испытаниях различных 
многоэлементных электрогенерирующих сборок [6]. 

В качестве иллюстрации возможностей создан- 
ного пакета программ рассмотрим эксперимен- 
тальные и расчетные характеристики пятиэлемент- 


ной профилированной ЭГС (длины ЭГЭ в ЭГС - 
9,55; 4,60; 4,00; 4,60; 9,55 см), основные параметры 
которой приведены в [7]. 


Таблица 3. Состав модулей для расчета тепловых процессов ЭГС 


Программный модуль 

Назначение модуля 


Расчет: 

псбі, пнеі, еш 

теплопроводности пара цезия в зазоре, 
гелия и изоляции системы теплосброса 

І5И 

теплового потока излучением 

сі5и 

распределения тепловыделения за счет 
нейтронного и гамма излучений 

КОЕК1, КОЕОЫ 

коэффициентов аппроксимации по ра- 
диусу и длине 

ТЕМРЕК 

температурного поля 

РАКСР 

средних параметров 


Таблица 4. Состав функциональных модулей для расчета ЭГС 


Функциональный 

модуль 

Назначение модуля 


Расчет: 

РЕС01 

температурного поля вакуумного и диффу- 
зионного режимов 

РЕС02 

вакуумных ВАХ 

РЕСОЗ 

переходного режима при заполнении ЭГС 
цезием 

РЕС04 

изотермических дуговых ВАХ 

РЕС05 

диффузионных ВАХ 

РЕС06, РЕС07 

статических и динамических ВАХ 

РЕС08 

переходного режима от диффузионного к 
дуговому с учетом поджига и гашения дуги 


Варианты распределения относительного те- 
пловыделения в топливных сердечниках ЭГС 
представлены на рис. 3. В соответствии с рекомен- 
дациями [6] для проведения расчетов выбран про- 
филь при <2ус=2. 



Рис. 3. Варианты распределения тепловыделения в пятиэле- 
ментной ЭГС 

На рис. 4 представлены теоретические кривые 
Рейзора для эмиттера (вакуумная работа выхода 
ГЦ= 5,0 эВ) и коллектора (І 7 С °=5,0 эВ), что соответ- 
ствует подложкам из вольфрама. ВАХ ТЭП (рис. 5), 
рассчитывались по эмпирической модели плазмы 
ТЭП [8] при межэлектродном зазоре ^=0,04 см, 
Р с = 266 Па, температуре коллектора Г с =900 К. 

На рис. 6 представлены расчетные (1, 4) и экс- 
периментальные (2, 5) кривые ВАХ пятиэлемент- 
ной профилированной ЭГС в штатном режиме ра- 
боты. 
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Рис. 4. Теоретические кривые Рейзора: 1) эмиттера и 2) кол- 
лектора 



Рис. 5. ВАХ ТЭП при варьировании температуры эмиттера Ту 
1) 1400; 2) 1600; 3) 1800; 4) 2000; 5) 2200; 6) 2300 К 



Рис. 6. ВАХ ЭГС в исходном состоянии и после изменения 
параметров модели. Тепловая мощность 0: 1-3) 2,5; 
4-6) 4,0 кВт 

В качестве иллюстрации дополнительных воз- 
можностей алгоритмического и программного 
обеспечения было предположено, что в ЭГС одно- 
временно произошли следующие изменения. В 
первом ЭГЭ уменьшилась теплопроводность кол- 
лекторной изоляции, во втором возникла утечка 
тока между электродами, в третьем снизились те- 
плопроводящие свойства изоляции между несущей 
трубкой и корпусом, а также возникла утечка тока 
через коллекторную изоляцию. 

Зависимости тока от напряжения нагрузки при 
данных предположениях представлены кривыми 3, 
6 на рис. 6. Основным достоинством разработан- 


ных программных модулей является не только рас- 
чет интегральных ВАХ ЭГС, но и распределений 
температур, потенциалов и токов в расчетных сече- 
ниях, которые позволяют провести анализ вну- 
треннего состояния ЭГС. В качестве примера на 
рис. 7, 8 представлено распределение температур, 
токов и напряжений в ЭГС при тепловой мощно- 
сти 0=4,0 кВт и напряжении на нагрузке С н =3 В. 



Рис. 7. Температуры в разных точках ЭГС: Т ю , Тп ~ в топлив- 
ных сердечниках по оси и на расстоянии радиуса пе- 
ремычки от центра, Т Е , Т с ~ эмиттера и коллектора, Т н , 
Т т - несущей трубки и теплоносителя 

Анализ графиков показывает, что в первом ЭГЭ 
повышена температура коллектора из-за уменьше- 
ния теплопроводности коллекторной изоляции. В 
третьем элементе наблюдается рост температуры 
несущей трубки и коллектора в связи со снижени- 
ем теплопроводности изоляции между несущей 
трубкой и корпусом. Во втором элементе при нали- 
чии электропроводящих мостиков между электро- 
дами за счет увеличения тока (рис. 8) происходит 
снижение уровня температурного поля эмиттера и 
разности потенциалов между электродами \] с . 



Рис. 8. Распределение температур эмиттера (Т Е ), коллектора 
(Т с ), токов (І с ) и напряжений (ІІ С ) в ЭГС при измене- 
нии параметров 


Наконец, в связи с утечкой тока через коллектор- 
ную изоляцию в третьем элементе происходит рост 
тока, снижение напряжения между электродами и 
уровня температур эмиттеров в последних трех ЭГЭ. 

Следует отметить, что электрический пробой 
коллекторной изоляции при принятой коммута- 
ции элементов в ЭГС приводит к повышению тем- 
пературного поля в ЭГЭ, расположенных до места 
пробоя и снижению температур эмиттера в ЭГЭ, 
расположенных после места пробоя относительно 
исходного уровня. Так в данном случае величина 
тока, проходящего через нагрузку по первому токо- 
выводу равна 9,6 А, по коммутационной перемыч- 
ке второго элемента 239,3 А, а по последним трем 
элементам и второму токовыводу - 141 ,7 А. Данное 
обстоятельство позволяет сделать вывод о необхо- 
димости перевода ЭГС в область короткого замы- 
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кания с целью гарантированного исключения пе- 
регрева эмиттера первого ЭГЭ [9] . По результатам 
измерения напряжения и токов на токовыводах 
предложен способ определения ЭГЭ, в котором 
произошел пробой коллекторной изоляции [6]. 

Полученные результаты свидетельствуют об 
адекватном качественном поведении выходных ха- 
рактеристик используемой математической модели 
при изменении теплофизических параметров и 
внутреннего состояния ЭГЭ в ЭГС. 

Разработанный комплекс программ нашел при- 
менение для оценки внутренних неизмеряемых пара- 
метров ЭГС, а также отработки гипотез при аномаль- 
но высоких временах переходных процессов [6, 10], 
получении гистерезисных статических ВАХ во время 
проведения испытаний энергонапряженных ЭГС по 
программе РКК "Энергия" на реакторе ВВР-К [6]. 

Выводы 

1. Создан пакет прикладных программ для моде- 
лирования термоэмиссионных систем, включа- 
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